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ABSTRACT 


This  was  a  two-part  numerical  study  using  ANSYS  Fluent  to  develop  novel  heat 
sink  concepts  from  first  principles.  The  objective  of  this  research  was  to  highlight 
geometric  structures  that  incorporate  the  principles  of  the  stack  effect  to  improve  the  heat 
transfer  capability  of  a  heat  sink  under  natural  convection.  The  first  part  investigated  the 
heat  transfer/fluid  flow  characteristics  of  vertically  aligned  tubes.  The  gaps  between  tubes 
break  up  the  thermal  and  velocity  boundary  layers  and  the  moving  fluid  within  a  tube 
entrains  the  cooler  ambient  air  surrounding  the  gap,  thus  increasing  mass  flow  rate  and 
average  Nusselt  number  through  each  tube.  The  optimal  gap-to-length  ratio  varies 
depending  on  the  number  of  tubes  in  the  system.  The  second  part  built  upon  the  insight 
gained  to  develop  heat  sinks  to  compare  to  pin-fin  heat  sinks.  A  tube  system  heat  sink 
provides  a  significant  improvement  in  the  heat  transfer  capability  over  a  circular  pin-fin 
arrangement,  demonstrated  by  an  increase  in  both  the  overall  heat  transferred  and  average 
heat  transfer  coefficient.  The  principles  discussed  in  this  study  have  the  potential  to 
expand  the  capability  of  natural  convective  heat  transfer. 
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I. 


INTRODUCTION 


This  thesis  conducted  a  two-part,  detailed  numerical  study  that  analyzed 
geometric  structures  that  have  the  potential  to  increase  the  overall  heat  transfer  capability 
of  a  heat  sink  when  incorporated  into  the  design.  The  first  part  of  the  study  explored  the 
possibility  of  increasing  the  mass  flow  rate  and  overall  heat  transfer  capability  through 
vertical  tubes  by  utilizing  principles  of  the  stack  effect  (i.e.,  using  multiple  inlets  to  create 
a  suction  effect  to  draw  cooler  ambient  air  into  the  tube).  The  second  part  of  the  study 
used  the  insight  gained  to  develop  novel  heat  sinks  and  compare  them  to  standard  pin-fin 
heat  sinks.  Three-dimensional  computations  of  flow  and  heat  transfer  were  carried  out  for 
constant  temperature  heat  addition  for  single  vertical  tubes,  multiple  stacked  tubes,  tube 
system  heat  sinks  and  circular  pin-fin  heat  sinks. 

A.  NATURAL  CONVECTION 

Natural  convection  heat  transfer  has  been  studied  quite  extensively  due  to  its 
many  applications  and  high  reliability.  In  natural  convection,  fluid  motion  is  induced  due 
to  density  variations  caused  by  a  temperature  gradient  in  the  working  fluid  that  is 
subjected  to  a  body  force  (e.g.,  gravity).  Elenbaas  [1]  thoroughly  explored  the  free 
convective  heat  transfer  from  the  inner  surface  of  vertical  tubes  of  different  cross-sections 
such  as  equilateral  triangle,  square,  rectangle,  and  circle.  More  recently,  Totala  et  al.  [2] 
experimentally  verified  that  the  local  heat  transfer  coefficient  is  a  maximum  at  the 
entrance  to  a  vertical  tube  subjected  to  free  convective  heat  transfer  and  subsequently 
decreases  as  the  thermal  boundary  layer  grows.  Davis  and  Perona  [3]  concluded  that  the 
development  entrance  length  for  free  convection  flows  in  vertical  tubes  is  quite  large  and 
the  assumption  of  fully  developed  flow  is  invalid  for  many  cases.  Due  to  the  work  of 
Davis  and  Perona,  it  is  imperative  to  understand  the  entrance  region  to  fully  grasp  the 
heat  transfer  capability  of  vertical  tubes  subjected  to  natural  convection.  Sieder  and  Tate 

[4]  developed  a  correlation,  eqn.  (1.1),  for  the  average  Nusselt  number  and  the 

Graetz  {Gz)  number  for  the  combined  entry  length  for  internal  laminar  flow.  The 
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combined  entry  length  means  that  both  the  thermal  and  velocity  boundary  layers  are 
developing  at  the  same  time. 


Nujj 


\0.14 


(1.1) 


0.60  <  Pr  <  5 


0.0044  < 


f  \ 

iL 


<9.75 


Nujj  >  3.66 

In  addition  to  Sieder  and  Tate,  Incropera  [5]  discusses  a  correlation  developed  by 
Hausen,  eqn.  (1.2),  for  the  thermal  entry  length  for  moderate  Prandtl  (Pr)  numbers, 
which  can  also  be  used  for  a  combined  entry  length  with  high  Pr  fluids.  Gnielinski  [6] 

developed  a  Nuj^  vs  Gz  correlation,  eqn.  (1.3),  for  fully  developed  turbulent  flows.  All 

three  of  these  correlations  are  discussed  by  Incropera  et  al.  [5]. 

-  0.0668  (D/L)Re„Pr 

Nu„  =3.66  + - ^ 

l  +  0.04[(D/L)Re„Pr]'  2) 

—  (//8)(Re„-1000)Pr 

”  l  +  12.7(//8)''^(Pr'=-l) 

B.  HEAT  TRANSFER  CAPABILITY 

When  improving  the  convective  heat  transfer  capability  from  a  surface  there  are 
three  possible  avenues:  increase  the  temperature  difference,  increase  the  heat  transfer 
surface  area,  or  increase  the  heat  transfer  coefficient.  The  most  widely  used  method  is  to 
increase  the  surface  area  by  the  addition  of  fins.  Different  designs,  concentrations,  and 
orientations  of  these  fins  have  been  analyzed  in  numerous  experimental  and  numerical 
studies. 


The  alignment  of  pin-fins,  either  staggered  or  inline,  was  shown  to  affect  the  heat 
transfer  performance  in  natural  convection  with  the  conclusion  that  inline  pin-fins  yield  a 
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higher  heat  transfer  rate  than  staggered  [7].  The  studies  conducted  by  Sertkay  et  al.  and 
Huang  et  al.  both  indicate  that  depending  on  the  design,  pin-fins  in  the  horizontal  and 
upward  vertical  direction  can  be  comparable  while  the  downward  vertical  direction 
provides  the  worst  heat  transfer  characteristics  [8,9]. 

During  the  literature  review,  only  two  studies  were  uncovered  that  analyzed 
hollow  pin-fins.  Elshafei  [10]  compared  a  hollow  pin-fin  to  solid  pin-fins  subjected  to 
natural  convection  under  different  orientations  and  concluded  that  heat  transfer 
performance  for  hollow  pin-fins  was  better  than  that  of  solid  fins,  which  resulted  in  a 
lower  temperature  difference  between  the  base  plate  and  the  ambient  air  for  a  given  heat 
input.  Awad  [11],  on  the  other  hand,  compared  straight  fins,  solid  pin-fins,  hollow  pin- 
fins,  and  convergent-divergent  fins  and  concluded  hollow  pin-fins  have  the  lowest  heat 
transfer  performance,  based  on  having  the  lowest  average  heat  flux.  In  both  of  these 
studies,  the  hollow  pin-fins  had  one  entry  point  at  the  base  of  the  pin-fin.  At  first  glance, 
these  two  articles  appear  to  contradict  each  other;  however,  this  study  offers  an 
explanation  as  to  how  both  are  correct  based  on  their  respective  experimental  setup. 

C.  NATURAL  VENTILATION 

For  large-scale  applications,  natural  convection  is  utilized  to  ventilate  buildings. 
There  are  three  general  techniques  for  natural  ventilation  depicted  in  Figure  1.  Namely: 
1)  single-sided  ventilation,  2)  wind-pressure  driven  cross-ventilation,  and  3)  buoyancy 
pressure-driven  stack  (chimney)  ventilation  [12]. 
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Figure  1.  Single-Sided  Ventilation  (Left);  Wind-Pressure  Driven  Cross- 
Ventilation  (Middle);  Buoyancy  Pressure-Driven  Stack 
Ventilation  (Right) 


In  stack  ventilation  single  or  multiple  chimneys  are  designed  with  various  air 
inlets  to  remove  the  hotter  internal  air  creating  a  suction  effect  that  draws  in  cooler  fresh 
air.  There  are  many  factors  that  influence  the  performance  of  stack  ventilation,  which 
includes  the  size  and  location  of  the  air  inlets  and  openings.  Krzaczek  et  al.  concluded 
that  the  stack  effect  could  be  generated  with  a  temperature  difference  as  little  as  12 
degrees  Celsius  [12]. 
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II.  NUMERICAL  MODEL  AND  DATA  REDUCTION 


The  computational  fluid  dynamics  software  utilized  for  all  models  discussed  in 
this  study  was  ANSYS  FLUENT.  At  the  basic  level,  the  problem  considered  for 
numerical  modeling  was  a  vertical  tube  under  natural  convection  conditions  subjected  to 
constant  temperature  heat  addition.  Many  additional  geometries  were  created  with 
multiple  tubes  stacked  vertically  with  a  gap  between  each  tube;  Figure  2  represents  the 
basic  geometries  analyzed.  The  left  tube  model  depicts  a  single  tube  with  a  length-to- 
diameter  ratio  equal  to  five  and  the  right  graphic  depicts  a  multiple  stacked  tube  model 
with  a  length-to-diameter  ratio  equal  to  one  and  a  gap-to-length  ratio  equal  to  0.4.  In  the 
multiple  stacked  tube  models  the  opening  is  a  gap  between  the  heated  tube  segments. 


Opening 


Adiabatic  Walls 


Adiabatic  Walls 


Adiabatic  Walls 


^  '  Adiabatic  Walls 


Opening 


Figure  2.  Representation  of  Basic  Geometries  Created  for  Initial  Analysis 
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In  all  geometries,  only  the  fluid  domain  was  analyzed,  which  simplified  the  model 
by  not  using  a  fluid-solid  interface  and  neglected  conduction  through  the  heat  sink. 
Figure  3  depicts  the  heat  sink  designs  analyzed. 


Figure  3.  Heat  Sink  Designs  Analyzed 

A.  MODEL  ASSUMPTIONS 

For  the  single  tube  and  multiple  stacked  tube  geometries,  the  vertical  tubes  were 
assigned  to  be  smooth  walled,  with  a  no-slip  boundary  condition,  and  maintained  at  a 
constant  temperature.  On  top  of  and  below  each  tube  was  a  box  that  contained  open 
boundaries  with  a  zero  Pascal  relative  pressure  across  the  boundary  with  an  ambient 
temperature  below  the  tube  wall  temperature.  The  horizontal  sections  directly  connected 
to  the  tubes  were  assumed  adiabatic,  smooth  walls  with  a  no-slip  boundary  condition  as 
shown  in  Figure  2. 
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The  heat  sink  surfaces  were  assigned  to  be  smooth  walled,  with  a  no-slip 
boundary  condition,  and  maintained  at  a  constant  temperature.  The  open  boundaries  were 
the  four  walls  and  top  of  the  fluid  space  with  a  zero  Pascal  relative  pressure  across  the 
boundary  and  an  ambient  temperature  below  the  tube  wall  temperature. 


All  numerical  models  were  solved  using  the  Boussinesq  buoyancy  approximation. 
Additional  assumptions  made  in  the  numerical  model  include: 

•  Fluid  flow  is  steady  state,  laminar,  and  incompressible. 

•  Thermodynamic -physical  properties  of  the  fluid  are  independent  of 
temperature. 

•  Viscous  dissipation  and  pressure  stress  terms  are  negligible. 

•  Thermal  radiation  effects  are  neglected. 

B.  GOVERNING  EQUATIONS 


Based  on  the  previously  described  assumptions,  the  governing  equations  for 
steady  fluid  motion  are: 


Continuity  :  V  •  a  =  0 

Momentum  :  {u  •  Vm)  =  -Vp  -l-pV^M  -fp^g 


f  p  _p  ^ 

V  y 


Energy  :  p„c^  (v  •  mT)  = 


(2.1) 

(2.2) 

(2.3) 


Where  the  third  term  on  the  right  hand  side  of  eqn.  2.2  is  the  buoyancy  force  from  the 
Boussinesq  approximation. 


C.  NUMERICAL  PROCEDURE 

The  governing  equations  were  solved  using  the  finite  volume  method.  The 
discretized  equations  were  solved  using  a  pressure-based  solver  via  the  SIMPLE 
algorithm.  Information  on  the  pressure-based  solver  utilized  by  ANSYS  FLUENT  can  be 
found  in  the  ANSYS  Fluent  Theory  Guide  and  the  SIMPLE  algorithm  is  discussed  by 
Patankar  [13,14].  For  the  single  tube  and  multiple  stacked  tube  geometries,  the  mesh  size 
was  refined  until  there  was  no  significant  change  in  the  results  from  one  mesh  size  to  the 
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more  refined  size.  The  residuals  were  driven  to  less  than  lE-6  and  convergence  was 
assumed  when  the  residuals  showed  negligible  change  over  100  iterations.  For  the  heat 
sinks,  a  fine  mesh  size  was  used  in  the  ANSYS  meshing  tool  and  convergence  was 
assumed  when  the  energy  residual  showed  negligible  change  over  100  iterations.  Further 
refinement  of  this  mesh  size  was  not  conducted  based  on  the  observations  while 
analyzing  single  tube  and  multiple  stacked  tube  geometries.  From  that  analysis  it  was 
determined  that  the  energy  residuals  will  converge  with  a  less  refined  mesh  than  the 
velocity  residuals.  It  was  noted  that  flow  patterns  did  not  change  with  a  refinement  of  the 
mesh.  The  change  was  in  the  velocity  of  the  fluid,  which  was  accompanied  with  an 
equivalent  change  in  the  temperature  difference,  resulting  in  a  negligible  change  in  the 
overall  heat  transfer.  Various  ANSYS  mesh  and  physics  reports  are  included  in 
Appendix  A. 

Since  the  only  heat  addition  into  the  system  was  via  the  heated  tube  walls,  the 
overall  energy  balance  was  calculated  as  follows: 


Energy  Balance = 


Energy^^,- Energy.^ 
Energy.^ 


*100  = 


(4^)  “(44L 

V  /  open  V  /  tube 


V  /  tube 


*100 


(2.4) 


The  average  energy  balance  for  the  single  tube  and  multiple  stacked  tubes  was  0.0027% 
and  the  energy  balance  in  the  worst  case  was  0.011%.  While  the  average  energy  balance 
for  the  heat  sinks  was  0.11%  and  1.72%  in  the  worst  case. 


A  heat  transfer  check  was  calculated  for  the  single  tube  and  multiple  stacked  tube 
models  by  comparing  the  heat  transfer  calculated  from  the  total  heat  generated  to  the 
summation  of  the  heat  transfer  from  each  individual  tube  as  follows: 


Heat  Transfer  Check  (%) 


100 


(2.5) 


The  two  methods  of  heat  transfer  averaged  a  variation  of  2. 12%. 


The  numerical  model  was  validated  first  by  solving  the  unsteady  laminar  flow 
condition  from  a  thermal  plume  such  as  a  cigarette,  additionally  the  model  was  used  to 
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solve  the  flow  and  temperature  field  for  a  tube  without  the  application  of  a  body  force 
(gravity)  to  ensure  that  no  fluid  motion  was  induced  and  subsequently  the  model  was 
solved  for  a  body  force  100  times  gravity  to  check  the  other  extreme  for  any  unexpected 
results.  With  the  numerical  model  satisfactorily  validated,  the  model  was  used  to 
calculate  the  flow  and  temperature  fields  for  vertical  tubes  of  various  length-to-diameter 

(L/D)  ratios.  The  calculated  was  plotted  against  \[g^  and  compared  to  the 
correlations  developed  by  Sieder,  Tate,  and  Gnielinski  [4,6]. 


D.  DATA  REDUCTION 


Both  dimensional  and  non-dimensional  analysis  methods  were  used  to  compare 
the  results  from  the  different  models.  The  average  heat  transfer  coefficient  of  each  tube 
was  calculated  as  follows: 


Qconv  = 


(2.6) 


(2.7) 


where 


(T  -T  ) 

y  m,i  m,o  J 


In 


^  'Y _ Y  ^ 

^  m.o 


T  -T 

\  s  mj  J 


(2.8) 


The  Nusselt  number,  ,  is  the  ratio  of  convective  to  conductive  heat  transfer  and  is 


defined  as: 


(2.9) 


thus 


NUr 


(2.10) 
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The  Prandtl  number,  Pr ,  is  the  ratio  of  kinematic  viscosity  to  thermal  diffusivity  and 
assumed  to  be  constant.  The  Reynolds  number,  Re  is  defined  as  the  ratio  of  inertial 
forces  to  viscous  forces  acting  on  a  fluid. 


RSd 


puD 

p 


(2.11) 


where 

=  (2.12) 


Therefore  the  Re^  can  be  re-written  as: 


Re^ 


4m 

nDp 


(2.13) 


The  Graetz  number,  Gz ,  is  used  to  determine  the  developing  lengths  of  flow  in  a  duct  and 
defines  as: 


Gz  = 


Pr 


(2.14) 


The  local  skin  coefficient,  Cj,  is  defined  as: 


C 


(2.15) 


The  Stanton  number,  St ,  is  a  modified  Nusselt  number  defined  as: 


St  = 


h 


puCp 


(2.16) 


The  numerical  models  neglected  radiation  even  though  radiation  has  the  potential 
to  provide  a  significant  portion  of  the  total  heat  transfer  when  using  natural  convection. 
Therefore,  the  radiative  heat  transfer  was  analytically  solved  and  included  for  analysis 
purposes  where  noted.  Figure  4  depicts  an  arbitrary  enclosure  used  to  determine  the 
effective  emissivity  of  a  partial  enclosure.  The  radiation  heat  transfer  was  calculated  as 
follows: 
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Arbitrary  Opening:  Area  =  A1 


Figure  4.  Arbitrary  Shape  and  Opening  Used  in  Radiative  Analytical  Analysis 


1 


1-1 


1-1 


V'^2  y  V 


A 


(2.17) 


(2.18) 


Therefore,  the  total  heat  transferred  was  the  summation  of  O  and  O  , : 

’  ^conv  ^racl 

Qlot  Qconv  Qmd  . 

The  thermal  resistance,  ,  for  the  evaluated  heat  sinks  was  determined  as: 

T  -T 

n  _  s  oo 
^Th  “ 


a 


(2.19) 


(2.20) 


The  thermal  entrance  length  is  the  distance  from  the  entrance  where  the  heat  flow 
is  fully  developed  (e.g.,  for  constant  temperature  heating)  Am  =3.66  and  is  calculated 
using  the  equations  discussed  by  Rathore  and  Kapuno  [15],  where  the  Re^  is  obtained 
from  eqn.  2.13. 

Laminar  Flow:  1 «  0.05  Re^  Pr  (2.21) 


Turbulent  Flow: 


(2.22) 
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The  thermal  entrance  length  was  then  used  to  calculate  the  percentage  of  the  thermal 
entrance  length  the  tube  utilized  for  flow  development  before  being  broken  up  by  the  next 
gap.  The  percentage  of  the  thermal  entrance  length  was  calculated  as  follows: 


%  of  Thermal  Entrance  Length  = 


—  *100. 

4 


(2.23) 


The  velocity  entrance  length  was  not  calculated  since  the  working  fluid  is  air  (Pr  = 
0.707)  the  thermal  boundary  layer  will  develop  much  faster.  The  thermal  entrance  length 
will  be  approximately  60%  that  of  the  velocity  entrance  length. 
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III.  RESULTS  AND  DISCUSSION 


A.  SINGLE  TUBE  AND  MULTIPLE  STACKED  TUBES  COMPUTATIONAL 
RESULTS 

The  numerical  computations  were  performed  for  a  range  of  governing  parameters 
shown  in  Table  1.  The  length-to-diameter(L/ D) ,  gap-to-length  (G/L)  and  number  of 

tubes  were  varied  and  all  additional  variations  are  a  direct  result  of  the  manipulation  of 
these  three  parameters. 


Tablet.  Governing  Parameters 


LID 

0.133-300 

GjL 

0.033-2 

Nujj 

6.1-24.3 

Re^ 

358-3151 

Gz 

7.4-4502.9 

Pr 

0.707 

T  -T 

s  00 

73  ‘’K 

I.  Effect  of  Length-to-Diameter  Ratio,  L/D 

For  a  single  vertical  tube,  the  length  of  the  tube  plays  a  vital  role  in  the  fluid  flow 
and  heat  transfer  characteristics  of  the  tube.  As  the  LID  ratio  increased,  the  mass  flow 
rate  through  the  tube  increases,  which  in  turn  increases  the  Re^.  As  the  Re ^ increases, 
the  thermal  entrance  length  for  the  tube  increases;  however,  as  shown  in  Figure  5,  the 
length  of  the  tube  increases  faster  than  the  thermal  entrance  length  therefore  a  larger 

percentage  of  the  entrance  length  was  used  by  the  flow,  which  resulted  in  a  smaller  Nuj^ 
for  the  tube. 
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Figure  5.  Percentage  of  Thermal  Entrance  Length  used  by  the  Developing 

Flow  vs.  Reynolds  Number 


When  plotting  the  determined  Nuj^  versus  yjcz  *  the  single  tube  models  trend 

with  the  Sieder  and  Tate  correlation  for  laminar  flow  in  a  combined  entrance  length  until 
the  flow  transitions  to  turbulent  [4].  Only  one  geometry  depicted  flow  that  is  fully 

turbulent  and  developed,  but  in  that  one  case  the  versus  Vg^  agrees  with  the 
Gnielinski  correlation  [6].  These  results  are  plotted  in  Figure  6. 
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Adapted  from  correlations  developed  by  Sieder  and  Tate,  Gnielinski,  and  Hausen  [4-6]. 


Figure  6.  Average  Nusselt  Number  vs.  Square  Root  Inverse  Graetz  Number 


2.  Effect  of  Stacked  Tubes 

Based  on  the  results  from  single  vertical  tubes  an  L  /  D  =  1  was  selected  for  the 
majority  of  the  investigations  into  stacking  multiple  tubes.  When  the  tubes  are  aligned 
vertically  the  moving  fluid  exiting  a  tube  is  at  a  slightly  lower  pressure  than  the  ambient 
air  surrounding  the  gap  causing  it  to  be  drawn  in  and  entrained  in  the  moving  fluid.  This 
forces  the  moving  fluid  from  the  previous  tube  to  utilize  a  smaller  cross-sectional  area  as 
the  entrained  ambient  air  develops  the  new  boundary  layer  thus  resulting  in  an  increase  of 
the  mass  flow  rate  through  each  additional  tube.  The  left  graphic  in  Figure  7  depicts  the 
fluid  streamlines  compressed  by  the  ambient  air  entrained  in  the  vertical  tube  system.  The 
right  graphic  depicts  the  same  streamline  with  the  addition  of  the  wall  heat  flux. 
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Figure  7.  Streamlines  of  Stacked  Tube  Model  with  Wall  Heat  Flux 


Interesting  trends  in  the  heat  transfer  characteristics  are  revealed  when  the  NUj^  is 

calculated  for  each  tube.  The  Nuj^  increases  through  each  tube  as  depicted  in  Figure  6 

plotted  alongside  the  single  tube  results,  the  Sieder  and  Tate  correlation  (combined  entry 
region),  the  Hausen  correlation  (thermal  entry  region),  and  the  Gnielinski  correlation 
(fully  developed  turbulent  flow)  [4-6] .  There  are  multiple  contributing  factors  that  result 

in  an  increase  of  the  Nuj^  . 
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•  The  thermal  and  velocity  boundary  layers  are  disrupted  due  to  the  gap 
between  tubes  therefore  the  boundary  layers  have  to  re-develop.  This  is 

shown  by  the  calculated  Nuj^  for  each  tube  trending  with  the  Sieder  and 

Tate  [4]  correlation  for  a  combined  entry  region.  This  can  also  be  seen  in 
Figure  7  by  the  streamlines  of  the  ambient  fluid  and  the  tube  wall  heat  flux 
contours. 

•  The  entrainment  of  the  surrounding  air  near  the  gap  increases  the  mass 
flow  rate,  which  in  turn  increases  the  length  of  the  thermal  development 
region.  Since  the  tube  length  is  maintained  constant  the  percentage  of  the 
entrance  length  utilized  decreases,  as  seen  in  Figure  5,  resulting  in  an 
increase  of  the  Nuj^ . 


The  local  skin  friction  coefficient,  C,  is  a  maximum  at  the  entrance  of  a 

round  tube  and  decreases  to  a  horizontal  asymptote.  The  skin  coefficient 
can  be  used  to  relate  the  wall  shear  stress  to  heat  transfer  through  the 
Cf 

Reynolds  Analogy,  — ^  =  St ,  which  can  be  used  to  solve  for  the  heat 


T  c  (T  -T^) 

transfer  coefficient  to  rewrite  the  heat  flux  as  ^ .  This 

u 


results  in  both  the  shear  stress  and  the  heat  transfer  being  a  maximum  at 
the  entrance  of  each  tube,  as  depicted  in  Figure  8. 
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Figure  8.  Stacked  Tube  Model  with  Associated  Wall  Shear  Stress  Plot 


3.  Effect  of  Gap-to-Length  Ratio,  G/L 

A  series  of  geometries  consisting  of  4,  8,  12,  16,  and  20-stacked  tubes  with  the 
Gj L  ratio  varied  from  1/30  to  2  was  analyzed.  One  additional  model  consisting  of  50- 
stacked  tubes  was  analyzed  to  assess  for  a  potential  maximum.  The  GjL  ratio  has  a 
significant  impact  on  the  heat  transfer  characteristics  of  the  system.  When  the  GjL  ratio 
is  too  small,  the  stacked  tubes  act  similar  to  a  single  tube,  where  the  Nuj^  through  each 
subsequent  tube  decreases.  That  trend  was  only  observed  when  the  G/L  =  1/30.  As  the 
Gj Lratio  is  increased  closer  to  2,  there  is  almost  no  degradation  in  the  amount  of  heat 
transferred  from  one  to  the  next. 

The  multiple  stacked  tubes  with  a  constant  LI D  =  I  diverged  from  the  Sieder 
and  Tate  [4]  correlation.  The  break  from  the  previously  established  trend  was  assessed  to 
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be  dependent  on  the  Re Figure  9  plots  Nuj^  versus  Remand  shows  an  inflection  point 
at  Re^  equal  to  1150. 
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Figure  9.  Average  Nusselt  Number  vs.  Reynolds  Number 


The  nature  of  the  shift  is  dependent  on  the  G/L  ratio.  When  the  G/L  =  1/3  the 
trend  appears  to  shift  down  and  follow  the  Hausen  [5]  correlation  for  the  thermal  entrance 
length  of  a  pre-developed  velocity  profile.  Figure  10  shows  the  velocity  profile  for  the 
first,  second,  ninetieth  and  twentieth  tubes.  It  is  observed  that  the  entrance  velocity 
profile  is  undeveloped  in  the  earlier  tubes  but  in  the  later  tubes,  the  entrance  velocity 
profile  shows  the  characteristics  of  a  developed  flow. 
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Figure  10.  Velocity  Profiles  for  20-Stacked  Tubes 

However,  for  a  G/L  =  2/3  the  divergence  breaks  in  a  different  direction.  Figure 
11  shows  the  entrance  and  exit  velocity  profile  for  the  fiftieth  tube,  which  shows  an 
atypical  entrance  velocity  profile.  The  entrance  velocity  profile  has  a  star  pattern  to  the 
velocity  distribution  but  by  the  exit  of  the  tube,  the  velocity  profile  displays  a  more 
typical  developed  velocity  profile.  This  star  pattern  appears  to  be  a  result  of  both  the  G/L 
ratio  and  the  tube  entrance  geometry. 
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The  additional  G/ L  ratios  evaluated  diverge  at  the  same  Re^ ;  however,  they  fan 
out  depending  on  the  G/L  ratio.  Figure  12  depicts  the  divergence  for  multiple  G/L 
ratios  but  these  results  were  obtained  from  models  with  a  mesh  not  as  refined  as  depicted 
in  Figures  6  and  9.  Sample  mesh  sizes  are  tabulated  in  Appendix  A. 
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Figure  12.  Average  Nusselt  Number  vs.  Reynolds  Number  for 

Less  Refined  Mesh 


The  optimal  G/L  ratio  was  determined  by  dividing  the  total  heat  transferred  from 
the  system  by  the  summation  of  the  length  of  the  individual  tubes  and  gaps  and  plotted 
against  the  G/L  ratio  for  the  different  stacked  tube  systems.  If  the  goal  is  to  maximize 
heat  transfer  in  a  small  overall  length  the  use  of  smaller  gaps  will  be  optimal;  however, 
for  larger  overall  lengths  the  use  of  larger  gaps  are  optimal  since  the  larger  gap  reduces 
the  degradation  of  heat  transfer  from  one  tube  to  the  next.  Therefore,  the  optimal  G/L 
ratio  is  dependent  on  the  overall  length  of  the  tube  system.  Figure  13  depicts  this  trend  by 
plotting  the  normalized  heat  transfer  per  unit  length  for  the  various  G/L  ratios. 


22 


1.2 


ti) 

c 

<D 

_l 

•ti  0.8 

c 

3 

.!U 


c 

(D 

i_ 

H 

-t-' 

ns 

Q) 

X 

"D 

OJ 


0.6 


0.4 

ro 

E 


0.2 


^.^4  Stacked  Tubes 
^►8  Stacked  Tubes 
Hirl2  Stacked  Tubes 
"^16  Stacked  Tubes 

H“20  Stacked  Tubes 


0 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8 

Gap-to-Length  (G/L)  Ratio 

Figure  13.  Normalized  Heat  Transfer  per  Unit  Length  vs.  Gap-to-Length  Ratio 
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The  results  of  the  heat  transferred  from  the  stacked  tubes  was  plotted  to  compare 
to  the  heat  transferred  from  a  single  tube  with  the  same  total  length.  From  Figure  14,  it 
becomes  apparent  that  if  the  GjL  ratio  is  too  large  the  heat  transfer  will  be  worse  than 
that  of  a  single  tube  with  the  same  total  length.  This  was  only  observed  in  the  4-stacked 
tube  systems,  since  the  GjL  ratio  where  the  overall  performance  becomes  worse 
increases  as  the  number  of  tubes  in  the  stacked  tube  system  increases.  It  is  also  observed 
that  a  stacked  tube  system  can  be  used  to  transfer  the  same  amount  of  heat  energy  with 
the  use  of  a  portion  of  the  length  that  would  be  required  of  a  single  tube,  in  some  cases 
the  stacked  tubes  would  only  need  -55%  the  total  length  of  a  single  tube  to  transfer  the 
same  amount  of  heat.  In  Figure  14,  the  heat  transfer  was  plotted  against  the  total  length  of 
the  system.  For  instance,  the  8-stacked  tubes  with  the  smallest  G/L  ratio  will  result  in  the 
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smallest  total  length;  therefore,  will  be  the  farthest  left  data  point.  This  is  the  same  for  all 
multiple  stacked  tube  systems  plotted. 


Figure  14.  Heat  Transfer  versus  Total  Length  of  the  System 


B.  TUBE  SYSTEM  VERSUS  PIN-FIN  HEAT  SINK  COMPUTATIONAL 

RESULTS 

Based  on  the  insight  gained  from  analyzing  stacked  tube  systems,  three  heat  sinks 
were  created  with  a  base  of  15cm  x  15cm  with  a  height  of  20cm.  Circular  pin-fin  heat 
sinks  were  created  with  the  same  spacing  and  outer  diameter  for  comparison.  Two 
additional  64  pin-fin  heat  sinks  were  created  to  have  a  surface  area  range  that 
encompassed  the  range  of  the  tube  system  heat  sinks.  Table  2  shows  the  surface  area  for 
each  heat  sink  evaluated  and  in  the  subsequent  figures  the  data  will  be  plotted  against  the 
surface  area. 
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Table  2.  Surface  Area  for  Each  Evaluated  Heat  Sink 


#  Eins 

1 

5 

9 

16 

25 

64 

64 

(mod) 

Pin-Ein  Surface  Area 
(m^) 

0.096 

0.154 

0.193 

0.299 

0.344 

0.382 

0.502 

Tube  System  Surface 
Area  (m  ) 

0.187 

0.303 

0.396 

0.476 

0.565 

N/A 

N/A 

1.  Effects  of  Pin-Fin  versus  Tube  System 

The  convective  and  radiative  heat  transfer  was  plotted  for  each  heat  sink  geometry 
versus  the  surface  area  in  Eigure  15.  Based  on  the  polynomial  curve  fit  the  9-  and  16-tube 
systems  show  an  18%  and  19%  improvement  in  the  convective  heat  transfer  over  an 
equivalent  surface  area  pin-fin  heat  sink.  The  9-tube  system  heat  sink  shows  a  118% 
improvement  over  the  9-pin-fin  heat  sink.  As  the  number  of  tube  systems  increases  the 
heat  transfer  improvement  diminishes  when  compared  to  an  equivalent  surface  area  pin- 
fin  heat  sink.  This  indicates  that  the  tube  system  geometry  has  an  optimal  point  and  loses 
its  effectiveness  as  the  number  of  tube  systems  continues  to  increase.  This  offers  an 
explanation  to  why  Aw  ad,  whose  experiments  consisted  of  95  staggered  or  100  in-line 
hollow  pin-fins,  concluded  that  hollow  pin-fins  perform  more  poorly  than  solid  pin-fins, 
while  Elshafei  used  8  hollow  pin-fins  and  determined  that  hollow  pin-fins  outperform 
solid  pin-fins  [10,11]. 
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Figure  15.  Heat  Transferred  vs.  Surface  Area 

Figure  16  depicts  the  average  convective  heat  transfer  coefficient  for  each  heat 
sink  and  clearly  shows  that  the  upward  oriented  tube  system  has  a  higher  average  heat 
transfer  coefficient  when  compared  to  pin-fins.  This  verifies  that  the  geometric  structure 
of  a  stacked  tube  system  enhances  the  overall  heat  transfer  characteristics  of  the  heat  sink 
and  not  only  due  to  the  additional  surface  area.  The  increased  heat  transfer  coefficient  is 
caused  by  two  effects:  1)  the  disruption  of  the  thermal  and  velocity  boundary  layers  due 
to  the  gaps  in  the  tubes  and  2)  the  creation  of  a  stack  effect  by  vertically  aligning  the 
tubes  that  causes  cooler  ambient  air  to  be  drawn  into  the  heat  sink 
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Figure  16.  Convective  Heat  Transfer  Coefficient  vs.  Surface  Area 


2.  Effects  of  Orientation 

The  various  heat  sinks  were  analyzed  in  three  orientations:  the  fins  pointing 
upward,  the  fins  pointing  downward,  and  the  fins  pointing  horizontal.  Figure  17  depicts 
the  convective  heat  transfer  versus  surface  area  for  each  of  the  orientations.  Based  on  the 
changes  of  the  heat  transferred  due  to  orientation  it  becomes  apparent  that  the  tube 
systems  are  more  greatly  affected  by  orientation  than  the  pin-fin  heat  sinks. 
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Figure  17.  Convective  Heat  Transfer  vs.  Surface  Area 


The  results  of  the  pin-fin  heat  sinks  correlates  with  results  previously  published 
by  Sertkaya  et  ah,  Huang  et  ah,  and  Elshafei  that  horizontal  pin-fins  can  be  comparable  to 
upward  oriented  pin-fins,  depending  on  the  pin-fin  arrangement,  while  the  downward 
facing  pin-fin  heat  sink  yields  the  lowest  heat  transfer  capability  [8-10].  Contrary  to  the 
pin-fin  heat  sinks,  the  downward  facing  tube  systems  provide  a  better  heat  transfer 
capability  than  the  horizontal  tube  systems.  Figure  18  depicts  the  normalized  heat  transfer 
based  on  orientation. 
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Figure  18.  Normalized  Heat  Transfer  Based  on  Orientation 


Figure  19  plots  the  thermal  resistance  of  the  different  heat  sinks  evaluated  at 
different  orientations  versus  the  surface  area.  This  shows  that  the  tube  system  heat  sinks 
have  a  lower  thermal  resistance  than  the  pin- fin  heat  sinks  and  the  16-tube  system  has  the 
lowest  thermal  resistance. 
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Figure  19.  Thermal  Resistance  vs.  Surface  Area  for  the  Various  Orientations 


3.  Effects  of  Additional  Appendages 

Based  on  the  streamlines  in  the  16-tube  system  heat  sink,  as  shown  in  the  upper 
left  and  right  of  Figure  20,  additional  features  were  evaluated  for  any  additional  heat 
transfer  enhancements.  The  placement  of  the  additional  appendages  is  shown  in  the  lower 
picture  of  Figure  20.  The  purpose  of  this  was  to  evaluate  if  impingement  cooling  could  be 
utilized  by  placing  an  extrusion  perpendicular  to  the  streamline  in  certain  locations.  Two 
additions  were  evaluated  both  separately  and  combined:  impingement  fins  and  a 
center  X. 
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Figure  20.  Fluid  Streamlines  and  Location  of  Additional  Appendages 

The  result  from  the  different  combinations  of  the  additions  reveals  a  reduction  in 
the  heat  transfer  capability  of  the  heat  sink  with  the  exception  of  the  downward  facing 
tube  system.  The  results  are  graphically  depicted  in  Figure  21. 
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Figure  21.  Comparison  of  Thermal  Resistances  for  the  Different  Additional 

Appendage  Combinations 
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IV.  CONCLUSIONS 


Heat  transfer  and  fluid  flow  characteristics  from  vertically  aligned  tubes  were 
examined  numerically.  The  effects  of  the  length-to-diameter  ratio,  number  of  stacked 
tubes,  and  gap-to-length  ratio  on  the  heat  transfer  and  fluid  flow  were  studied. 

The  numerical  results  closely  follow  the  previously  established  correlations  of 
Sieder  and  Tate  for  the  combined  entry  length  under  laminar  flow  conditions.  It  was 
observed  that  the  gap  between  the  tubes  breaks  up  the  thermal  and  velocity  boundary 
layers.  The  cooler  ambient  air  near  the  gap  becomes  entrained  in  the  fluid  exiting  the  tube 
causing  the  mass  flow  rate  and  average  Nusselt  number  to  increase  through  each 
subsequent  tube.  The  effects  on  the  heat  transfer  and  fluid  flow  characteristics  are 
dependent  on  the  number  of  stacked  tubes,  gap-to-length  ratio,  and  the  length-to-diameter 
ratio.  It  was  numerically  demonstrated  that  a  stacked  tube  system  with  a  total  length  55% 
of  a  single  tube  could  transfer  the  same  amount  of  heat. 

The  insight  gained  was  utilized  to  develop  novel  heat  sinks  to  compare  to  circular 
pin-fin  heat  sinks.  The  effect  of  number  of  fins  and  orientation  on  the  heat  sinks  was 
additionally  analyzed.  The  tube  system  heat  sink  outperforms  the  pin-fin  heat  sink; 
however,  that  is  dependent  on  the  number  of  tube  systems  utilized.  If  the  number  of  tube 
systems  is  too  great,  the  performance  will  actually  be  worse  than  an  equivalent  surface 
area  pin-fin  heat  sink.  The  16-tube  system  heat  sink  provided  the  lowest  thermal 
resistance  and  highest  average  heat  transfer  coefficient  out  of  all  of  the  heat  sinks 
analyzed.  The  orientation  of  a  tube  system  heat  sink  is  more  drastically  affected  by  the 
orientation  than  pin-fins  but  contrary  to  the  pin-fin  heat  sinks  the  downward  facing  tube 
system  outperforms  the  horizontal  tube  system  with  the  upward  facing  tube  system 
offering  the  best  heat  transfer  capability.  An  upward  facing  tube  system  heat  sink  has 
been  shown  to  provide  a  20%  increase  in  the  heat  transfer  capability  over  an  equivalent 
area  pin-fin  heat  sink. 

The  results  of  this  numerical  study  show  that  there  are  geometric  structures  that 
can  be  incorporated  to  improve  the  overall  heat  transfer  capability  of  a  system  subjected 
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to  natural  convection.  In  this  specific  study,  a  tube  system  heat  sink  was  shown  to 
outperform  a  circular  pin-fin  heat  sink.  There  is  still  further  progress  that  can  be  made  by 
optimizing  the  spacing  between  tube  systems,  the  alignment  of  the  systems  (e.g.,  in-line 
versus  staggered),  and  the  entrance  and  exit  geometry  of  each  tube  in  the  tube  system. 
An  experimental  validation  of  these  results  was  conducted  with  the  results  discussed  in 
Appendix  B. 

Upon  the  conclusion  of  this  study,  Plunkett  Associates  used  a  Direct  Metal  Laser 
Sintering  process  to  construct  novel  heat  sinks  that  could  not  be  created  using  traditional 
manufacturing  processes  [16].  These  novel  heat  sinks  demonstrated  a  significant 
improvement  in  heat  transfer  over  typical  extruded  heat  sinks.  The  importance  of  this 
press  release  validates  the  assumption  that  unused  geometric  structures  exist  that  can  be 
utilized  to  increase  the  natural  convection  capability  of  heat  sinks  and  companies  are 
starting  to  explore  the  possibilities. 
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APPENDIX  A.  ANSYS  MESH  AND  PHYSICS  REPORTS  EOR 

VARIOUS  GEOMETRIES. 

A.  ANSYS  PHYSICS  REPORT 


Table  3.  20-Stacked  Tubes  ANSYS  Physics  Report 


Domain 

Boundaries 

Boundary  -  open_boundary 

Type 

PRESSURE-OUTLET 

chimney  20tubes 

Boundary  -  tube  (heat  source) 

Type 

WALL 

Boundary  -  wall  chimney 20tubes  (adiabatic  walls) 

Type 

WALL 

B.  ANSYS  MESH  REPORT 


Table  4.  Various  ANSYS  Mesh  Reports 


Domain 

Nodes 

Elements 

20-Stacked  Tubes 

942,870 

5,213,424 

50-Stacked  Tubes 

2,616,515 

14,569,406 

16-Stacked  Tubes  (less  refined) 

149,497 

742,468 

9  Tube  System  Heat  Sink 

416,397 

2,391,599 

64  Pin-Lin  Heat  Sink 

385,639 

1,930,144 
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APPENDIX  B.  EXPERIMENTAL  VERIEICATION. 


To  assess  the  validity  of  the  numerical  results  two  heat  sinks  were  created  using 
some  of  the  principles  uncovered  during  the  numerical  analysis.  The  heat  sinks  were 
designed  using  a  square  pin-fin  for  ease  of  manufacturability.  A  0.25”  end  mill  bit  was 
used  to  cut  the  slots  between  the  square  pin-fins,  or  the  hollow  square  pin-fin  heat  sink  a 
drill  press  with  a  3/16”  drill  bit  was  used  to  hollow  out  the  pin-fin  and  additionally  to 
provide  two  entrance  regions  into  the  hollowed  out  fin.  The  manufacturing  of  the  heat 
sinks  were  lined  up  via  eyeball  vice  using  a  computer  controlled  script  and  the  vice  had  a 
little  bit  of  play  in  the  track  system  that  held  it  to  the  table.  Figure  22  depicts  the  heat 
sinks  created  for  verification. 


Figure  22.  Heat  Sinks  Manufactured  for  Experimental  Verification 


The  heat  sinks  were  modeled  in  SOLIDWORKS  accounting  for  imperfections 
created  during  the  machining  process.  The  heat  sinks  were  assessed  using  the  same 
numerical  model  used  in  previous  analysis. 
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A.  EXPERIMENTAL  SETUP 

For  the  experimental  verification,  a  medical  grade  1”  by  1”  WATLOW  ceramic 
heater  was  selected  as  the  heat  source.  The  heater  has  one  type  K  thermocouple  to 
measure  the  temperature  at  the  interface  between  the  heater  and  the  heat  sink.  Each  heat 
sink  was  evaluated  at  four  power  settings  ranging  from  0.11  to  1.54  Watts  and  the 
temperature  was  recorded  at  each  power  setting  once  the  temperature  stabilized.  The 
experiment  was  conducted  in  a  lab  space  with  an  ambient  air  temperature  of  23.0  C  at  the 
time  of  the  experiment. 

The  numerical  models  were  analyzed  with  both  a  constant  temperature  and 
constant  heat  flux  wall  boundary  condition  on  the  heat  sink.  The  trends  in  both  cases 
followed  the  experimental  results.  The  constant  heat  flux  cases  were  utilized  for  analysis 
since  the  results  more  closely  followed  the  magnitude  of  the  experimental  data. 

B.  DISCUSSION 

Figure  23  depicts  the  surface  temperature  contour  with  an  ambient  fluid 
Streamline  of  the  hollow  square  pin-fin  heat  sink. 
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Figure  23.  Hollow  Square  Pin-Fin  Heat  Sink  Surface  Temperature  with  Air 

Streamline  Overlay 


Figure  24  is  a  comparison  of  the  heat  flux  versus  temperature  difference  for  the 
experimental  and  numerical  results.  This  shows  that  the  hollow  square  pin-fin  heat  sink 
has  a  lower  average  heat  flux  for  a  given  temperature  difference  when  compared  to  the 
solid  square  pin-fin  heat  sink.  This  result  correlates  with  Awad’s  [11]  conclusion  that 
hollow  pin-fins  have  the  lowest  heat  transfer  performance,  based  on  having  the  lowest 
average  heat  flux. 
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Temperature  Difference  (C) 

Figure  24.  Comparison  of  Heat  Flux  vs.  Temperature  Difference 

Figure  25  is  plot  of  the  temperature  difference  versus  the  power  input  from  the 
heater.  This  figure  shows  that  the  hollow  square  pin-fin  heat  sink  has  a  lower  temperature 
difference  for  a  given  power  input.  This  result  correlates  with  Elshafei’s  [10]  findings 
that  hollow  pin-fins  have  better  heat  transfer  performance  than  solid  pin-fins  resulting 
from  a  lower  temperature  difference. 
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Figure  25.  Comparison  of  Temperature  Difference  vs.  Power  Input 

Based  on  the  previous  numerical  results  in  this  study  one  possibility  to  reckon  the 
difference  between  Aw  ad  [11]  and  Elshafei  [10]  was  based  on  the  optimal  number  of 
hollow  pin-fins  in  the  heat  sink,  but  Figures  24  and  25  also  show  that  depending  on  how 
their  data  was  plotted  they  could  have  determined  different  conclusions. 

The  average  convective  heat  transfer  coefficient  is  compared  to  the  temperature 
difference  in  Figure  26.  This  figure  shows  that  both  experimentally  and  numerically  the 
solid  square  pin-fin  heat  sink  will  have  a  higher  heat  transfer  coefficient  that  the  hollow 
square  pin-fin  heat  sink.  This  result  was  not  expected  based  on  the  previous  numerical 
analysis  conducted  with  circular  pin-fin  heat  sinks. 
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Figure  26.  Comparison  of  Convective  Heat  Transfer  Coefficient  vs. 

Temperature  Difference 


-♦-Square  (Experimental) 
-■-Square  (ANSYS) 

Hollow  (Experimental) 
— Hollow  (ANSYS) 


The  hollow  square  pin-fin  heat  sink  does  show  a  lower  temperature  difference  for 
the  same  power  input  when  compared  to  the  solid  square  design,  which  correlates  to  a 
smaller  thermal  resistance  as  shown  in  Figure  27.  The  lower  thermal  resistance  is  due  to  a 
larger  heat  transfer  surface  area  vice  an  increase  in  the  heat  transfer  coefficient. 
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Figure  27.  Comparison  of  Thermal  Resistance  vs.  Temperature  Difference 


There  are  a  few  significant  conclusions  that  can  be  drawn  from  this  experimental 
validation. 

•  The  experimental  and  numerical  data  trend  together.  Therefore 
experimentally  validating  the  numerical  model  used  in  the  previous 
analysis  of  stacked  tube  system  heat  sinks. 

•  Depending  on  how  data  is  plotted  different  conclusions  can  be  made. 
Therefore,  it  is  vital  to  look  at  the  results  from  multiple  angles  and 
assimilate  all  the  data  to  make  a  proper  conclusion. 

•  The  hollow  square  pin-fin  heat  sink  has  a  lower  thermal  resistance  and 
uses  less  material  than  the  solid  heat  sink.  Therefore,  the  hollow  pin-fin 
has  better  overall  heat  transfer  characteristics  than  a  solid  pin-fin,  despite 
having  a  lower  average  heat  transfer  coefficient. 

•  The  cross-sectional  profile  of  the  fin  will  have  a  significant  impact  on  the 
performance  of  the  hollow  pin-fin.  The  numerical  study  with  circular 
hollow  pin-fins  shows  an  increase  in  the  average  heat  transfer  coefficient 
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in  addition  to  increasing  the  heat  transfer  surface  are  whereas  the  square 
hollow  pin-fin,  both  numerically  and  experimentally,  shows  a  decrease  in 
the  average  heat  transfer  coefficient  when  compared  to  its  solid  pin-fin 
counterpart. 
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